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ABSTRACT 

Analytic distribution functions (dfs) for the Galactic disc are discussed. The DFs de- 
pend on action variables and their predictions for observable quantities are explored 
under the assumption that the motion perpendicular to the Galactic plane is adiabat- 
ically invariant during motion within the plane. A promising family of dfs is defined 
that has several adjustable parameters. A standard df is identified by adjusting these 
parameters to optimise fits to the stellar density in the column above the Sun, and 
to the velocity distribution of nearby stars and stars ~ 1 kpc above the Sun. The 
optimum parameters imply a radial structure for the disc which is consistent with 
photometric studies of the Milky Way and similar galaxies, and that 20 per cent of 
the disc's luminosity comes from thick disc. The fits suggest that the value of the V 
component of the Sun's peculiar velocity should be revised upwards from 5.2 kms 
to ~ 11 kms . It is argued that the standard df provides a significantly more reliable 
way to divide solar-neighbourhood stars into members of the thin and thick discs than 
is currently used. The standard df provides predictions for surveys of stars observed 
at any distance from the Sun. It is anticipated that dfs of the type discussed here will 
provide useful starting points for much more sophisticated chemo-dynamical models 
of the Milky Way. 



Key words: 

hood 



galaxies: kinematics and dynamics - The Galaxy: disc - solar neighbour- 



ly 1 1 INTRODUCTION 



A major thread of current research is work directed at 
understanding the origin of galaxies. There are excellent 
prospects of achieving this goal by combining endeavours in 
three distinct areas: observations of galaxy formation taking 
place at high redshift, numerical simulations of the gravita- 
tional aggregation of dark matter and baryons, and stud- 
ies of the Milky Way. The latter field is dominated by a 
series of major observational programs that started fifteen 
years ago with ESA's Hipparcos mission, whic h returned 
paral laxes and proper motions for ~ 10 s stars (|Perrvmanl 
1 19971 ) . Hipparcos established a more secure astrometric ref- 
erenc e frame, and the UCAC2 catalogue (Zac harias et al.l 
|2004 ) uses this frame to give proper motions for sev- 
eral million stars. These enhancements of our astrometric 
database have been matched by the release of m ajor pho- 
tometric catalogues [ DENIS (Epc htein et alj|2005h . 2MASS 
l|Skrutskie et alj|2006l h SDSS (|Abazaiianll2009h ] and the ac- 
cumulation of enormous numbers of stell ar spectra, start- 
ing with the Geneva-Copen hagen Survey (|Nordstrom et ail 
12004 iHolmberg et al.ll2007l . hereafter GCS) and continuing 
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with the SPSS, SEGU E (|Yannv et all l2009bh and RAVE 
l|Steinmetz et alj[2006h surveys - on completion the SEGUE 
and RAVE surveys will provide 0.25 x 10 6 and ~ 10 6 low- 
dispersion spectra, respectively. These spectra yield good 
radial velocities and estimates of [Fe/H] and [a/Fe] with er- 
rors of ~ 0.1 dex. Two surveys (HERMES and APOGEE) 
are currently being prepared that will obtain large numbers 
of medium-dispersion spectra from which abundances of sig- 
nificant numbers of elements can be determined. The era of 
great Galactic surveys will culminate in ESA's Gaia mis- 
sion, which is scheduled for launch in late 2011 and aims to 
return photometric and astrometric data for 10 9 stars and 
low-dispersion spectra for > 10 7 stars. 

The Galaxy is an inherently complex object, and the 
task of interpreting observations is made yet more difficult 
by our location within it. Consequently, the ambitious goals 
that the community has set itself, of mapping the Galaxy's 
dark-matter content and unravelling how it was assembled, 
can probably only be attained by mapping observational 
data onto sophisticated models. We are developing a mod- 
elling strategy that has as its point of departure analytic ap- 
proximations to the distribution functions (dfs) of various 
components of the Galaxy (McMillan et al. in preparation) . 
In this paper we present such approximations for the thin 
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and thick discs. The paper is organised as follows. Section 
2 explains how the df is assembled. Section 3 compares the 
DF's predictions for various observables to data. In particu- 
lar evidence is presented that the Sun's V velocity is con- 
ventionally underestimated by ~ Gkms^ 1 and predictions 
are given for velocity distributions as a function of distance 
from the plane. Evidence is presented that the standard df 
provides a cleaner division of solar-neighbourhood stars into 
members of the thin and thick discs than has been available 
hitherto. Section 4 sums up and looks ahead. 



2 CHOICE OF THE DF 

Our approach to Galaxy modelling starts from dfs that are 
analytic functions of the action integrals (J r , J$, J^>) of or- 
bits in an integrable, axisymmetric Hamiltonian (McMillan 
et al., in preparation). The action associated with the az- 
imuthal invariance of the Hamiltonian is simply the z com- 
ponent of angular momentum L z , and the action denoted Jo 
in Binney & Tremaine (2008; hereafter BT08), which quan- 
tifies motion perpendicular to the symmetry plane z — 0, 
is here conveniently denoted J z . We shall be largely con- 
cerned with orbits that have sufficiently large values of L z 
that a reasonable approximation to their dynamics can be 
obtained by considering the motion parallel to the plane to 
proceed regardless of the vertical motion, and the vertical 
motion to be affected by the motion in the plane only in 
as much as the latter causes the force perpendicular to the 
plane to vary in time slowly enough for the vertical motion 
to be adiabatically invariant - see BT08 §3. 6. 2(b) for a jus- 
tification of this approximation. At any radius R we define 
the vertical potential 



$ z {z) = ${R,z)-${R,0), 



(1) 



where <&(R, z) is the full potential. Motion in <!> z has the 
energy invariant 

E„(z,v,) = \vl+$ z {z). (2) 

Given a value for E z , the vertical action can be obtained 
from a one- dimensional integral 

2 



J Z (E Z ) 



dz v z 



(3) 



where $ z (z m ax) = E z . 

Motion parallel to the Galactic plane is assumed to be 
governed by the radial potential 



= *(.R,0), 
so the radial action is 
1 f R * 

J r (E R ,L z ) = - / drv R , 
71 Jr„ 



(4) 



(•>) 



where R p and R a are the peri- and apo-centric radii and 
vr = y/2(E R - - Li/R?. 

Given a point in phase space, we can evaluate L z = J$ % 
E z and Er and thus obtain J r and J z , so a given df can be 
evaluated at any point in phase space. 

We start from the simplest plausible dfs, which have 
the form 




Holmberg Flynn 00 
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Figure 1. Comparison of the predictions of two theoretical mod- 
els and the density of A (triangles) and F (squares) stars versus 
distance from the plane from Holmberg & Flynn (2000). The full 
curve is obtained from the DF with cr z = 6.3kms — 1 while the 
dotted curve is the classical exponential p oc exp(— <E>/(j^) for the 
same velocity dispersion. The gravitational potential is that of 
Model II of 82.7 in BT08. 



Here /i is primarily responsible for determining the surface 
density of the disc, f r controls the degree of epicyclic motion 
within the disc, and f z controls the disc's vertical structure. 
Since there is a close relation between a star's angular mo- 
mentum L z and the radii to which it contributes to observ- 
ables, the appearance of L z in f r and f z makes it possible 
for the disc to become hotter and/or thinner at small radii. 



2.1 Vertical profiles 

We now consider the form of the function f z in equation (|6]) , 
which controls the disc's vertical structure. We focus on mo- 
tion in the solar cylinder of stars for which L z ~ Rqv c (Ro) 
so they do not make large radial excursions. For these stars 
f z is effectively a function of only J z . The classical choice 
of df is that of an isothermal sheet f z ~ e~ Ez ^ az l|Spitzerl 
1 19421 ). Our modelling strategy requires that we eliminate 
E z in favour of J z . In a separable potential Q, Z J Z = (i> z ), 
where fl z = 8E z /dJ z is the vertical frequency and the angle 
brackets denote a time average along the orbit with action 
J z . Moreover, by the virial theorem (vf) = E z in a har- 
monic oscillator, so we replace E z with Q Z J Z and arrive at 
what we shall refer to as the "pseudo-isothermal" df 



/a.(J«) 



27r/ oo d.7.e- I W' 



(7) 



where the denominator ensures that f ax satisfies the nor- 
malisation condition 



dz dv z f az — 1 ■<=> 



= to- 



rn 



f(Jr,J Z ,L z ) = fl(L z )f r {Jr,L z )f z (J z ,L z ). 



(6) 



In general a z is a function of L z to control the scale height 
as a function of radius, but for the moment we neglect this 
dependence and investigate the vertical density profile pre- 
dicted by the df by taking § z {z) to be the potential 
above the Sun in Model II of §2.7 in BT08; this model is 
disc dominated. 

The full curve in Fig. [T] shows the density profile pre- 
dicted by the df (O for a z = 6.3 km s , while the dotted 
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Figure 2. Upper panel: the full curve shows the vertical density 
profile predicted in the potential above the Sun by the DF of 
equation J9} with 7 = 2.6 and Vj = 18.7kms • The dotted 
curve is the profile predicted by equation J5J with a z = 18 kms - 1 . 
The points show the density of main-sequence stars measured 
by Gilmore & Reid (1983), Kuijken & Gilmore (1989) and Juric 
et al. (2008). Lower panel: the full and dashed curves show the 
distributions in v z at z = predicted by the DFs of equations (J9j 
and J7J, while the histogram shows the distribution in v z of the 
GCS stars. 



curve shows the classical isothermal p oc exp( — &/a z ). The 
two curves are very similar because equation J7J predicts 

that (vf^ 1 2 moves in a narrow range from 6.55 kms -1 at 
z — to a peak value 6.7 kms -1 at z = 240 pc. Both pre- 
dictions are in reasonable agreement with the densit ies of A 
and F stars measured bv lHolmberg fc Flvnnl (|2000T > shown 
by triangles and squares, respectively. 

The dotted curve in the upper panel of Fig. [2] shows 
the vertical density profile predicted by equation when 
o z = 18 kms -1 . The curvature of the profile has the wrong 
sign to fit the density profile o f dwa rfs with 4 < My < 
5 measured by iGilmore fc Reidl (|l983f ) , which is shown by 
circles. The velocity dispersion in this model is {t> z } 1 '' 2 = 
19.2 ± 0.3 kms -1 independent of z, so the DF is very close 
to an isothermal. The dotted curve in the lower panel shows 
that the Gaussian distribution in v z at z = predicted by 
the DF is a poor fit to the distribution in v z of the GCS 
stars. 

One way to obtain a vertical profile that is steeper at 
small heights and flatter at large heights is to replace the 
exponential in equation (0 with an algebraic function of 
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Figure 3. Velocity dispersion as a function of distance from the 
plane as predicted by the DF l|9)l . 



J z . The full curves in Fig. [2] show the density profile and 
velocity distribution predicted by the DF 



2tt J °° d j z (n. 



(9) 



with 7 = 2.6 and V~, = 18 . 7 km s . The agreement with the 
data of IGilmore fc Reidl <|l983T) is essentially perfect. This 
fit, and all subsequent fits, were obtained by adjusting the 
parameters by hand and judging the quality of the fit by 
eye. 

The curve in Fig. [3] shows the extent of the increase in 
velocity dispersion with height that is required to produce a 
thick-disc like flattening in the density profile at z > 500 pc: 
(^z) 1 ^ 2 rises from 19.6 kms -1 at z = to 37.5kms -1 at 
z — lkpc and 50.4 kms -1 at z — 2kpc. Also shown are 
three sets o f data points: filled squares show the values of a z 
inferred bv lKuiiken fc Gilmore! (11989T ) (hereafter KG89) for 
the velocity dispersion of K dwarfs; filled triang les show the 
disper sion of metal-poor M dwarfs inferred by iFuchs et al.l 
(2009); open pentagons show the analyt ic fit to the disper- 
sions of disc stars that was published bv lBond" et al.l t 2009). 
Sadly, the data do not tell a coherent story. Each data set 
is for a different stellar population, so in principle their dif- 
ferent trends in a z (z) could be matched by different density 
profiles. Therefore in Fig. [2]we plot the de nsity profiles asso - 
ciated with the Bond et al. sample (from ljuric et al1l2008l ) 
and the KG89 sample. We see that the Juric et al. density 
profile is less steep than the one from KG89, which is in- 
consist ent with the highe r velocity dispersions measured by 
KG89. IFuchs et al.l (|2009t ) do not give a density profile, but 
it would be surprising if the kinematics and vertical struc- 
ture of the M dwarfs were fundamentally different from that 
of the K dwarfs studied by KG89, so the extremely rapid 
increase in their values of a z (z) shown in Fig. [3] is hard to 
understand. 

The fact that the model curve in Fig.0agrees best with 
the KG89 may reflect the fact that the potential in which the 
DF is evaluated was constrained to be compatible with value 
for the surface densi ty of material that l ies wi thin 1.1 kpc 
of the plane given bv lKuiiken fc Gilmore! (|l99ll ), which was 
based on the KG89 data. In a gra vitat ional potent i al tail ored 
for the data of lJuric et all l|200SD and lBond et~ai1 (|2009l ) the 
DF might reproduce the data in these papers better than that 
of KG89. We do not pursue this question here. 
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Despite the simplicity of the df ((9J and the accuracy 
with which it fits the data, we will not employ it further 
because it does not provide a decomposition of the disc into 
populations of different ages, and there is no natural way 
of incorporating it into a df that also describes the radial 
structure of the disc. 

Star formation is known to have continued in the disc 
throughout the life of the Galaxy and the velocity dispersion 
of any cohort of coeval stars is known to increase secularly 
as a result of scattering by spiral arms and molecular cloud s 
|Spitzer fc Schwarzschildlll953l : ICarlberg fc Serrwoodlll985l) . 



Moreover, as the Galaxy ages, the chemical composition of 
the stars that are forming at a given radius changes, so stars 
formed at different times and different radii are in principle 
distinguishable. Hence it is useful to consider the disc's ag- 
gregate df to be a sum of the dfs for stars of different ages 
and velocity dispersions. 

We assume that the df of stars of age r is the "pseudo- 
isot hermal" df (|7|) with a 7 increasing with r according to 
(e.g. lAumer fc Binnevll2009l ) 



O"z0 ( — 



+ T1 



(10) 



Here <t z o is the velocity dispersion of stars at age r m ~ 
10 Gyr, n sets velocity dispersion at birth, and /3 — 0.38 
is an index that determines how a z grows with age. If we 
further assume that the rate of star formation has declined 
with time as e - *^* , then the aggregate df will be 



/thn(Jz) = 



j;-dre T /*o/ CTz (J,) 



(11) 



t ( e W*a _ 1) ' 

where a z depends on r through equation (|1U|) . In Fig. [3] 
we plot the vertical density profile produced by this aggre- 
;ate df when to = 8 G yr, ri = 0.1 Gyr, a z o = 20kms~ 1 
Aumer fc Binnevl [2009l l . We see that with these parame- 
ters we obtain a reasonable fit to the thin disc. Specifically, 
at 2 > 150 pc the density is nearly exponential with a scale 
height of 255 pc. In view of the dramatic difference between 
this pure thin-disc structure and the thin plus thick disc 
structure furnished by the algebraic df ([9]), it is perhaps 
surprising that the velocity distribution in the lower panel 
of Fig. [4] differs as little as it does from the dashed curve 
in the lower panel of Fig. [2] This comparison illustrates an 
important point: thick-disc stars spend relatively little time 
near 2 — so they contribute only inconspicuous wings to 
the velocity distribution there. The velocity distribution pre- 
dicted by /th n is close to Gaussian: the velocity dispersion 
rises from 16.5 kms -1 at 2 = to 19.6 kms -1 at 2 = 500pc 
and 20.9 kms -1 at z — lkpc. 

Fig. E] shows that a perfect fit to the Gilmore & Reid 
(1983) measurements can be obtained by adding a pseudo- 
isothermal component with a z = 38 kms -1 to the thin disc 
shown in Fig. [4] Within this structure (v"^ 1 ^ 2 increases from 
19 kms -1 at 2 = to 39 kms -1 at 2 = 1.5 kpc and then very 
slowly increases to 41.5 kms -1 at 2 = 2.5 kpc. It is worth 
noting that adding the thick disc increases the scale height 
in the exponential fit to the profile at low 2 from 255 pc to 
336 pc. 



2.2 Profiles within the plane 

[Shul (|l969T ) discussed dfs for planar discs of the form 
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Figure 4. Full curves: the vertical density profile (upper panel) 
and the distribution of v z (lower panel) predicted by the compos- 
ite DF mil l. Dashed line: an exponential of scale height 255 pc. 
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Figure 5. The result of adding to the composite thin disc of Fig. [4] 
a pseudo-isothermal disc with a z = 38 kms -1 that contains 20 
percent of the total mass. The dashed line is an exponential with 
scale height 336 pc. 



fr(E,L z ) 



-(B-E c )/ai 



(12) 



where E C (L Z ) is the energy of a circular orbit of angular 
momentum L z and a r (L z ) is a function that determines the 
velocity dispersion in the disc as a function of radius. By 
analogy with the vertical df we could rep l ace E — E c b y 
n r J r , where fi r = dE/dJ r (|Binnevl [l987l ; iDehnenl ll999T ) . 
However, the decrease in Q r as J r — -> 00 is so rapid that 
for sufficiently eccentric orbits the product il r J r decreases 
with increasing J r . Consequently, if one substitutes f2 r J r for 
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E — E c , at fixed L z and large J r the df increases with J r . 
To prevent this unphysical behaviour we replace E — E c by 
kJ t , where 



lim tt r (J r ,L z ) 



(13) 



is the epicycle frequency. Hence in this paper we adopt as 
the planar df of a pseudo-isothermal population 



fr{Jr, L z ) — e 



iJ r /o 



(14) 



In view of the normalisation condition |[8} , the mass dM 
placed by the df ((6| on orbits with angular momentum in 
the range (L z , L z + dL z ) is 



dM = (27r) 2 /id^ 



d J r e 



tJr/o 



= (2w) 2 ^fidL z (15) 



In the limit a r — > of a cold disc, only circular orbits are 
populated and this mass is equal to the mass 2-nY.R dR in the 
annulus (R,R + dR), where E is the disc's surface density. 
Hence for a cold disc 



fi(L*) = 



kRY, dR 
2nai dL z 



(16) 



where Rc(L z ) is the radius of the circular orbit of angu- 
lar momentum L z and the second equality uses the identity 
dL z /dR = Rn 2 /2n. Here we consider the case of an expo- 
nential disc, E = Y, e- (R - Ro)/Rd , where R d ~ 2.5 kpc is the 
scale length of the disc and Rq ~ 8 kpc is the Sun's distance 
from the Galactic centre. We assume that ay declines expo- 
nentially in radius with a scale length that is roughly twice 
that of the surface density 

o8 (iJ -Hc)/Kd 



o r {L z 



a r0 e" 



where q ~ 0.5. 



(17) 



This choice is motivated by naive epicycle theory, which im- 
plies that with q ~ 0.5 the scale height will be constant 
|van der Kruit fc Searle|[l982l ') provided a z /a r = constant. 

A df such as /i times equation (I14|l that is an even func- 
tion of L z does not endow the Galaxy with rotation. We in- 
troduce rotation by adding to the df an odd function of L z , 
which will not contribute to either the surface density or the 
radial velocity dispersion. A convenient choice for this odd 
contribution to the df is tanh(L z /Lo) times the even contri- 
bution, where Lo is a constant that determines the steepness 
of the rotation curve in the central region of solid-body ro- 
tation. At radii so large that Rv c 3> Lo this choice for the 
odd part of the df simply eliminates counter-rotating stars. 
We choose Lo = lOkms -1 kpc, a value sufficiently small for 
counter-rotating stars to be confined to the inner kiloparsec 
of the Galaxy, which is in reality bulge-dominated. Hence 
we consider the "pseudo-isothermal" planar df 



f(T r (Jr , L z 



ftE 



[1 +tanh(L z /L )]e' 



-nJ r /<r$ 



(18) 



It is interesting to evaluate the observables predicted 
by the df (1181) when the circular speed is a power law in R, 
v c = v (R/R)) a/2 . Then 



* R {R) = ^(R/Ro 



(19) 



In the limit a — > of a perfectly flat circular speed 
$ R {R) = vo\n{R/R ). (20) 




R/kpc 



2 4 6 8 10 

R/kpc 

Figure 6. From top to bottom, the surface density, rotation speed 
and radial velocity dispersion obtained from the DF l)18j l when 
the gravitational potential is given by equation d20D with vo = 
220kms — 1 and the velocity dispersion function is given by ||17|I 
with a>o = 28.6 km s — 1 and q = 0.45. In the top panel the dashed 
line shows an exponential with scale length 2.5 kpc; in equation 
(llSI t i?d = 2.25 kpc was used but the recovered surface density 
profile is clearly shallower than the exponential with this scale 
length, which is shown by the dotted line. 



Fig. [6] shows the surface density, rotation curve and ra- 
dial velocity-dispersion profile predicted by the pseudo- 
isothermal df ()18p for a flat circular-speed curve with vo = 
220 km s -1 . For the plotted profiles the function E(i? c ) has 
been taken to be an exponential of scale length 2.25 kpc, 
while the full curve in the top panel shows that the surface 
density produced by this choice of E(_R C ) provides a good 
approximation to the surface density of an exponential disc 
with a longer scale length, 2.5 kpc, which is shown by the 
dashed line. At the price of replacing the analytic function 



Efflc) 



3 (flo-Sc)/Hd 



(21) 



with a tabulated fu nction, the sur face density can be made 
exactly exponential (|Dehnenlll999h . Here we adopt the sim- 
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Figure 7. Comparison of the observed density of GCS stars in the 
( V R, V 4>) plane with the prediction of the planar pseudo-isothermal 
DF dl8l with parameters as in Fig. [6] The smooth curves show the 
projections of the model density onto the Vj, and vr axes, while 
the histograms show the corresponding distributions of stars in 
the GCS. 



pier expedient of using a slightly smaller value of Rd than 
the scale length of the disc we wish to produce. 

In the middle panel of Fig. [6] the mean rotation speed 
declines from 225 km s -1 at R = lkpc to ~ 200 km s -1 at 
R — 5kpc before slowly rising to 210kms _1 at 7? = lOkpc. 
The bottom panel shows that the radial velocity dispersion 
declines throughout the disc as expected, being ~ 34kms~ 1 
at Ro. 

Fig.[7]compares the distributions of v^> and vr velocities 
predicted by the pseudo-isothermal DF with the correspond- 
ing distributions of GCS stars with heliocentric velocities 
converted to vj, and vr assuming that the circular speed 
is 220kms _1 and the S un's velocity is vr = — lOkms" 1 , 
v,fi = v c + 5.2 km s -1 (|Dehnen fc Binnevl fl998. hereafter 
DB98). The theoretical and observed distributions in vr are 
satisfyingly similar, but the theoretical v$ distribution lies 
above the observed distribution at small v$ and well below 
it at Vtf, ~ 200kms _1 . We return to this issue below. 

As discussed in Section 12.11 realistically we must con- 
sider the thin disc to be a superposition of pseudo- isothermal 
cohorts of different ages and chemical compositions. The 
temperature of the pseudo-isothermal DF (|18[) is set by the 
parameter oyo through equation (|17[) . By analogy with equa- 
tion (fTTJ)) . we should make a r age dependent by adding to 
the right side of this equation the appropriate function of r. 
Then oy is given by 



a r (L z ,r) 



oyo 



a{R -R c )/R d 



(22) 



With o r (L z ,r) thus defined, it is natural to consider the 
thin-disc DF 



t ( e W*o - l) 



(23) 



where f tTr is defined by equation (|18p with a r (L z ,r) now ob- 
tained from equation (122[) . The DF /( J r , L z ) correctly yields 
values of observables averaged through the thickness of the 
disc. In the case of an external galaxy such averaged observ- 
ables are of interest, but samples of Milky- Way stars rarely 
if ever provide a sample that is unbiased in z. In particular, 
stars with large J z spend little time near the Sun so samples 
of local stars are biased against them. Since a star with large 
J z is likely to be old, it is likely to have large J r also. Hence 
samples of local stars are biased towards small J r and the 
DF (|23[) cannot be used to predict the properties of solar- 
neighbourhood stars, or indeed stars at that lie within any 
restricted range in z. Instead we must use the full thin-disc 

DF. 



3 FULL DF 

Putting together the planar and vertical parts of the DF for 
the thin disc introduced above, we have 



j thn( Jr , J z, L z ) — 



i (e T »/*o - 1) 



(24) 



where f ar is defined by equations (|18jl and l|22p . and f {Jz is 
defined by equation (0 but with Q z and a z now functions 
of L z through R c . By analogy with equation (|22p we have 



c z (L z ,t ) = <J Z Q 



n 



a(Ro-R c )/R d 



(25) 



An orbit's vertical frequency £l z is a function of all three 
actions, J r , J z and L z . However, the Jeans theorem assures 
us that the DF remains a solution of the collisionless Boltz- 
mann equation if in Q. z we set J r = 0. Restricting the J r - 
dependence of the DF in this way makes the DF easier to work 
with and is consistent with the work of Section 12.11 which 
was restricted to orbits that have J r — ("shell orbits"). 
Therefore in the following we do this. 

Our final thin-disc DF (|24[) is characterised by eight free 
parameters, Lo, Rd, q, oyo, o"z0, fl, T\ and r m . 

3.1 Thick disc DF 



At the end of Section \2. II we saw that the observed vertical 
density profile at the Sun can be reproduced by adding to 
the composite DF of the thin disc a pseudo-isothermal com- 
ponent with a z o = 38kms _1 that contains 20 percent of the 
mass. This result suggests that we add to the DF (|24|) of the 
thin disc the thick-disc DF 



/thk( Jr, Jz,L z ) = f ar (J r , L z )f< Tz (Jz), 



(26) 



where f (7r and f (Tz are defined by equations (|18[1 and (O 
above with ay and a z given by 



cr z (L z 



' (Jzo e 



q(R -Rc)/R d 



(27) 
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Table 1. Parameters of the standard DF (upper section) and val- 
ues used by Bensby et al. (2003) (lower section) 





Thin disc 


Thick disc 


Lo 


lOkms 


10 kms 


u 


z.zo kpc 


2.3 kpc 




0.45 


0.45 


O" r-0 


33.5kms _1 


60 km s — 1 


O"z0 


19 kms -1 


32 kms" 1 






0.24 


P 


0.33 




to 


8Gyr 




Tl 


110 Myr 




T m 


lOGyr 




<?U 


35 kms -1 


67 km s — 1 


<jy 


20 kms" 1 


38 kms" 1 


cr w 


16 kms -1 


35 km s — 1 


Vk 


15 kms -1 


46 km s — 1 



» Fuchs + 09 
o Bond + 09 



80 - 



40 - 



20 



z/kpc 

Figure 8. The full curve shows the radial velocity dispersion 
in the standard model. The dotted and dashed lines show the 
dispersions in the thin- and thick-disc components, respectively. 
The triangles show estimates of or for metal-poor M dwarfs from 
Fuchs et al. (2009), while the open pentagons show the analytic 
fit to o-flO) given by Bond et al. (2009) 



and in equation (|18p we use (cf eq. I21|l 

thke ( 



s = fe ,, Q (-Ro-«o)/«d 



(28) 



Here fc t hk is the ratio of thick to thin-disc stars in the solar 
cylinder. In principle this thick-disc DF introduces a further 
six parameters: Lo, Rd, q, ovo, & z o and fcthk- Table[T]lists all 
the parameters of the standard DF. We have not explored the 
option of using a different value of Lo for each disc because 
this parameter has a negligible impact on comparisons with 
local data. 

The choice of <r z o and k t hk for the thick disc is straight- 
forwardly made by the r equirement that the vert ical density 
profile match the data of lGilmore fc Reidl (1 19851 ). The choice 
of Rd, q and oyn for the thick disc is more problematic. 
Clearly these parameters should be constrained by the ra- 
dial density and kinematics of the disc at z > 1 kpc, where 
the thick disc is domi nant. The s t ronge st constraints are 
provided by the SDSS. lJuric et all (120081 ) found the scale- 
length of the thick disc to be similar to that of the thin disc. 
iBond et all l|2009l ) give an analytic fit to the dependence of 



1 kpc 



model 

model+error 



100 

v ( /km s 

1.75 kpc 




Figure 9. Distributions of v$ velocities lkpc (top panel) and 
1.75 kpc (lower panel) above the plane. The data points are from 
Ivezic et al. (2008) for stars with 0.8 kpc < z < 1.2 kpc (upper 
panel) and 1.5 kpc < z < 2 kpc (lower panel). The full curves show 
the model distributions at z = lkpc and 1.75 kpc. The dashed 
curves show the result of convolving this distributions with the 
measurement errors of Ivezic et al. 



cjr on z out to z ~ 4 kpc, and iFuchs et al.l (l2009f) give sey - 
eral values of or at z < 800 pc. Finally Ivezic et all (|2008T 1 
provide distributions in v$ in several ranges of z. Figs [8] and 
show these data. 

Although the data sets are less clearly inconsistent than 
they are in Fig. [3] the data from Fuchs et al. clearly show 
a significantly steeper gradient than those from Bond et al. 
In the model or{z) has a slope intermediate between these 
values, and agrees with the data at z<l kpc. At greater 
heights it lies above the Bond et al. data, just as the model 
curve does in Fig. [3] 

The model's predictions for the distribution in v,/, at 
z = lkpc and z — 1. 75 kpc are shown in Fig. [9] together 
with da ta points from llvezic et al.l (|2008l ) the heliocentric 
data of llvezic et alj (|2008l ) have been converted to galac- 
tocentric velocities assuming ?^(Q) = 225.2kms _1 , which 
arises because in our adopted potential the circular speed 
is 220 kms " 1 and the peculiar V velocity of the Sun is 
5.2 kms" 1 dDehnen fc Binnevl[l998h . The full curves are the 
true model velocity distributions, and the dotted curves 
show the result of convolving these distributions with the 
errors reported by Ivezic et al., which are 19.5 kms -1 at 
z = lkpc and Sl^kms" 1 at z = 1.75 kpc. 

The model curves fall below the data at < 50 km s~ 
because in this region halo stars dominate the data points 
and the model is for the disc alone. Elsewhere the dotted 
curves provide a moderate fit to the data points. The fit 
in the upper panel would be improved by moving the data 
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points a few kms -1 to the right, which would be the effect of 
the upward revision of the Sun's peculiar velocity advocated 
below. The model curves are slightly too broad. Reducing 
the parameter oy in the thick-disc df makes them narrower, 
but this change also shifts the model curves still further 
to the right, and thus makes the overall fit less good. The 
distributions can also be made narrower by either increasing 
the thick-disc scalelength Rd or by decreasing q. However 
either change decreases the importance of stars at apocentre 
(which have angular momentum L z < L z (0)) relative to 
those at pericentre and thus exacerbates the predicted excess 
of stars at large v^. Extensive experimentation suggests that 
the fits shown in Fig. [9] cannot be significantly improved 
upon with a df of the form (|26[) . 

Because they are extracted from proper-motion data, 
the observational distributions in Fig. [9] are sensitive to the 
photometric distances employed. A possible resolution of the 
conflict in the lower panel of Fig. [5] between the model and 
data is that the distances employed are slightly too small: 
using larger distances would increase heliocentric velocities 
and thus cause the observational points to move away from 
the Sun's assumed velocity, = 225.2 kms" 1 . Another pos- 
sible resolution of the conflict between data and models in 
Fig. [5] is the increasing inaccuracy of the assumption of adi- 
abatic invariance of the vertical motion as random motions 
become more important. In a future publication this possi- 
bility will be examined with models based on orbital tori. 

3.2 The standard DF and the solar 
neighbourhood 

Fig. HOI shows prediction of the standard df for the structure 
of the solar neighbourhood. Full curves are for the whole 
disc and dashed curves show the contribution of the thin 
disc. The upper panels show for stars seen in the plane the 
distributions in v^, and vr after integrating over the other 
two velocity components. Comparison of these panels with 
the panels of Fig. [7J is instructive. The distribution of vr 
velocities is in reasonable agreement with the data, while 
the distribution of v^ velocities of Fig. [TO] agrees with the 
data better than the distribution in v$ in Fig. [7J Two fac- 
tors contribute to the improved fit to the velocities. First 
introducing a sum of pseudo-isothermals enhances both the 
core and the wings of the distribution - this effect is appar- 
ent in the distributions of vr velocities. More significantly, 
including vertical motion suppresses the distribution at 
low t>0 because this wing of the distribution is populated by 
stars with small values of R c that reach the solar neighbour- 
hood because they have large random velocities. On account 
of those large velocities, they have low probabilities of being 
found close enough to the plane to be included in the GCS. 
These stars are most likely to be observed near apocentre, 
when they have low values of vr, so depressing the contri- 
bution to the GCS of such stars does not suppress the wings 
of the distribution of vr velocities. 

The plots shown in Fig. [10] are obtained by setting the 
function E(i? c ) for the thin disc that appears in the df to 
the same exponential with scale length 2.25 kpc that was 
used to obtain Figs [6] and [Jj The use of 2.25 kpc in E(i? c ) is 
important not only to ensure that the disc's surface density 
is approximately exponential with the larger scale length 
2.5 kpc, but also to ensure that the predicted distribution of 



i>0 velocities agrees with the GCS data: when E has scale 
length 2.5 kpc, the predicted distribution falls off too slowly 

at v<t, > v c . 

The model v$ distribution deviates from the data in two 
respects: it lacks the pronounced peak in the data centred 
on Vfi = 205 kms -1 , and it extends too far on the high- 
velocity wing. The first shortcoming undoubtedly reflects 
the axisymmetry of the model and is discussed below. The 
second shortcoming , whic h is also evident in the fits to the 
data of llvezic et alj (|2008h for the thick disc (Fig. [9}, is more 
interesting. It can be moderated by increasing the param- 
eter q of equation (|22|) . This parameter controls the radial 
gradient of the thin disc's velocity dispersion, and a rapid 
decrease in the amplitude of epicycle motions at R c > Ro 
limits the number of stars with large angular momentum 
that can reach the Sun and thus depopulates the high-t;^ 
wing. However, an increase in q boosts the model profile at 
low so the overall agreement with the data is not im- 
proved unless the parameter oyo of equation (|22[) is simulta- 
neously decreased, and such a decrease leads to the model 
vr distribution being narrower than the data warrant. 

Oort's relation e.g., BT08 eq. (4.317) 



-B 



A-B 



(29) 



implies that the width of the model v^ distribution can be 
decreased relative to that of the vr distribution by changing 
from a flat to a falling circular-speed curve. However, one 
finds that the relative narrowing of the v^ distribution that 
is produced by adopting the power-law potential (|19|) with 
a = —0.2 produces a negligible improvement on the fit for 
constant circular speed shown in Fig. 1101 

The bottom left panel of Fig. 1 10 1 shows that the overall 
df p rovides an excellent fit to the vertical density profile 
from lGilmore fc Reidl (|l983l ). and that the vertical profile of 
the thin disc is extremely close to exponential. The latter 
result appears to be fortuitous in that it involves a subtle 
interplay between the non-trivial vertical force law and the 
number of stars with large random velocities that visit the 
solar neighbourhood from significantly nearer the Galactic 
centre. 

The bottom right panel of Fig. llOl shows that in both the 
thin and thick discs the mean rotation speed declines with 
distance from the plane. In the plane the thin disc rotates 
faster than the thick disc, as one naively expects. However, 
the rotation rate of the thin disc declines faster with z than 
that of the thick disc. The slower decline in the thick disc 
arises because we have set R& = 2.3 kpc in the function 
E(7? c ) = e~ Rc/Rd for the thick disc - with R d = 2.25 kpc for 
both discs the thick disc rotates ~ 15 kms -1 slower than the 
thin disc at all values of z. We have chosen i?d = 2.3 kpc for 
the thick disc to obtain a better fit to the long-dashed curve 
in the lower right panel of Fig. 1101 which is an analytic fit to 
the mean rotation rate extract e d from the proper motions 
of SDSS stars by llvezic et~ail (|200Sh . The dependence of 
the rotation rate on _R<j is consistent with the Stromberg 
equation 



2v c 



din 7? 



o% dz 



(30) 



For our preferred DF the asymmetric drift of the thick 
disc increases from only 20 kms -1 at z = to 90 kms - 
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Figure 10. Structure at the solar radius predicted by the standard DF (eqs, I24l and 1261 1. Full curves are for the entire disc while dashed 
curves show the contribution of the thin disc. The upper panels show velocity distributions for 2 = marginalised over the other velocity 
components. The vertical dotted line marks the local circular speed. In the lower panels dotted curves show the contributions of the 
thick disc. In the bottom-right panel the long-dashed line is the empirical fitting function of Ivezic et al. (2008). The values of the df's 
parameters are given in Table 1. 



at Z = 3kpc, consistent with th e values 30 — 50kms 
usually reported by observers (e.g. lEdvardsson et all 1 19931 : 
iGilmore et allll989l ). 

Fig. [TT] compares the observed density of stars in the 
plane of velocities (U, V) with respect to the LSB0 (top 
panel) to that predicted by the standard DF (bottom). The 
dynamic range in density that can be sampled with the 
GCS stars is limited, so only a portion of the standard 
df's predictions for the (U, V) plane is tested. Moreover, 
the limited number of GCS stars leads to the steepness of 
density gradients being underestimated, for example around 
(U, V) = (O^Ojkms- 1 . 

Near (U,V) = (0,0), the observational diagram shows 
density enhancements, or "streams", that are not bounded 
by the roughly ellipsoidal surfaces in velocity space on which 
actions are constant. Consequently, by Jeans' theorem, the 
presence of streams indicates that either the Galaxy's po- 
tential is not axisymmetric, or the Galaxy is not in a steady 
state - no DF that is a function of actions only can repro- 
duce these streams, although one hopes to be able to re- 



1 We follow BT08 (p. 12) in defining the LSR to be v R = v z = 
and = v c (Ro). The LSR is sometimes taken to be the speed 
of a closed orbit through Rq . In the presence of ephemeral spiral 
structure this second definition is probably not useful. 



produce them by perturbing such a DF using Hamiltonian 
perturbation theory. The streams account for much of the 
disagreement between the theoretical and observed veloc- 
ity distributions in Fig. [TO] In particular they account for 
the peak in the observed distribution lying ~ 15kms 
below the circular speed. 

3.3 The solar motion 

We have seen that the DF has difficulty simultaneously fit- 
ting the observed distributions in Vr and v^> of the GCS 
stars (top panels of Fig. I10p . A related problem was en- 
countered in the fit to the proper motions of thick-disc stars 
(Fig. [9}. The agreement between theory and data in both 
figures would be improved by shifting the observed v$ dis- 
tribution to the right. Such a shift would correspond to in- 
creasing the Sun's peculiar velocity by a few kilometers per 
second. 

Although such a shift significantly exceeds the formal 
error of 0.6 km s -1 on V© given by DB98, we should not 
lightly dismiss the possibility that Vq has been underesti- 
mated. The value given by DB98 was obtained by extrapo- 
lating to zero velocity dispersion a plot of V velocity versus 
squared velocity dispersion S 2 for stars grouped by colour 
such as that shown in Fig. 1121 Stromberg's equation (|30[1 
suggests that this relation will be linear if the square bracket 
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Figure 11. The density of GCS stars in the (U,V) plane (up- 
per panel) and that predicted by the standard DF (lower panel). 
Colours indicate log 10 of the stellar density and the horizontal 
line is at the circular speed. The upper plot was obtained by ap- 
plying the FiEstas algorithm of Ascasibar & Binncy (2005) to 
the GCS data and setting the density to zero if there was no star 
within 25kms _1 of a point. 



is constant, and Fig. 1121 shows that the Hipparcos data are 
consistent with this expectation if the groups with the low- 
est velocity dispersions (shown in red) are discounted. How- 
ever, Stromberg's equation is obtained under the assump- 
tion that the Galactic potential is axisymmetric, so in the 
limit of vanishing velocity dispersion, stars move on circu- 
lar orbits. In reality the potential has a non-axisymmetric 
component of amplitude ~ 7kms _1 , which manifests it- 
self, inter alia, by causing a plot of terminal veloc ity versus 
Gala c tic longitude to undulate at this level (e.g. iMalhotral 
ll995l ; lBinnev fc Merrifieidlll998l , Fig. 9.16). Given the non- 
axisymmetric component of the potential, linio-^o {V) can 
differ from the velocity of circular motion by of order the 
amplit ude ~ 7kms~ 1 of the no n-axisymmetric component. 

As lolling fc Dehnenl l|2003l ) pointed out in their deter- 
mination of the Oort constants, the larger the velocity dis- 
persion a population has, the less it is likely to be affected 
by non-axisymmetric forces associated with spiral struc- 
ture. The non-axisymmetric potential of the Galaxy's bar is 
thought to be responsible for the "Hercules stream" , an over- 
density of stars in velocity space around = v c — 50 kms , 
but there is no evidence that it significantly perturbs the 
velocity distribution at v^>v c , where the standard model 
conflicts with the data. Consequently, Fig. [10] offers an op- 




„ 1000 „ 
S /(km s ) 

Figure 12. The data points show, for each group of 1000 main- 
sequence stars ordered by colour, the solar motion versus the 
group's squared velocity dispersion in the plane of the sky. The 
points are based on the proper motions obtained by van Leeuwen 
(2007) from a reanalysis of the Hipparcos telemetry for samples of 
stars defined by Aumer & Binney (2009). The straight line is the 
least-squares fit to the black data points, which has y intercept at 
5.05 km s — 1 . The black rectangle on the y axis shows the proposed 
solar motion Vq = llkms -1 . The dashed line shows the solar 
motions predicted by model DFs for populations ranging in age 
from 1 to 11 Gyr under the assumption that Vq =11 kms" 1 and 
S 2 ~ (1.2a fl ) 2 (DB98). 
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Figure 13. Upper panel: the cumulative distributions of stars 
with smaller than the given value for the GCS data (jagged 
curve) that the standard DF under the assumption Vq = 
5.2kms _1 from DB98. Lower panel: the full curve shows the 
Kolmogorov— Smirnov probability that the distribution of V ve- 
locities of GCS stars is drawn from a model distribution as a 
function of the amount AV added to 5. 2 kms -1 . Dotted curve: 
the same but after randomly redistributing those observed stars 
that have velocities with respect to the conventional LSR in the 
range (-20, 10) kms" 1 . 
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portunity to determine Vq that is at least as valid as the 
traditional route using Stromberg's equation. 

The jagged curve in the upper panel of Fig. [13] shows 
the cumulative distribution in of the GCS stars under the 
assumption that v c (Ro) = 220 kms" 1 and Vq = 5.2 kms" 1 . 
The smooth curve shows the cumulative distribution of the 
model shown in Fig. 1101 The need to shift the distribution 
of GCS stars to the right is evident. 

The full curve in the lower panel of Fig. [13] shows the 
Kolmogorov-Smirnov probability Pks that the distribution 
of V velocities of GCS stars is drawn from the model shown 
in Fig. [10] as a function of the amount AV added to the 
solar motion given in DB98. For all choices of AV, Pks is 
small, largely due to the impact of streams on the data for 
— 20 km s _1 < Vlsr ^ 10 kms -1 . The impact of streams on 
Pks can be reduced by randomly redistributing stars within 
this range of Vlsr- The dotted curve in Fig. 1131 shows the 
dependence of Pks on AV when the randomised sample is 
compared to the model. The peak in Pks rises by three 
orders of magnitude and shifts from AV = 6.9 kms -1 to 
5.8 kms" 1 . 

Thus the stars that should be least affected by spiral 
structure suggest that DB98 underestimated Vq by AV = 
5.8 kms -1 , so the true solar motion is Vq = 11 kms -1 . The 
systematic error ~ lkms" 1 on this value is clearly much 
greater than the formal error reported by DB98. In Fig. [T2] 
the black rectangle marks this revised solar motion. 

It is interesting to test the extent to which Stromberg's 
equation is verified by pseudo-isothermal dfs for main- 
sequence stars of a given colour. The dashed curve in Fig. 1121 
shows the relation between S 2 ~ (1.2<tr) 2 (DB98) and V a 
that one obtains by calculating these quantities for a df of 
the form (|24[) with r m increasing from 1 to 11 Gyr; over this 
age range an increases from 12.4 to 28 kms" 1 . The dashed 
curve is plotted on the assumption that the solar motion is 
Vq = 11 kms" 1 , as marked by the black rectangle. At low 
S 2 the slope of the dashed curve is similar to the slope of 
the observational relation, but the slope flattens perceptibly 
with increasing S 2 . This flattening implies that the square 
bracket in Stromberg equation (|30[) diminishes with increas- 
ing velocity dispersion. The is no evident reason why this 
bracket should be constant. 

It is not inconceivable that spiral structure and/or the 
bar have shifted the observational points with S 2 in the 
range (250, 700) downwards from a relation that runs from 
the black square, between the red points and on to just above 
the points at S 2 > 1200. Moreover, a proponent of the con- 
ventional value of Vq should worry that if the dashed curve 
were moved down to start at that value of Vq, it would lie 
below nearly all the data points. We conclude that although 
we cannot confidently recommend an upward revision of Vq, 
considerable caution should be exercised in the use of the 
conventional value and more work is needed on the effect 
that spiral structure has on the local velocity space. 

Analysis of the space velocities of 18 ma ser sources 
for w hich trigonometric parallaxes are available ( Rcid ct al. 
2009) provides tent ative support for Vq bein g revised up- 
wards to 11 kms" 1 (|McMillan fc Binnevl[2009l ). 

The valu es of Vq, the proper mot ion of Sgr A*, 6.38 ± 
0.04 mass" 1 l|Reid fc Bru nthalcr 2004), a nd the distance to 
Sgr, A*, 8.33 ±0.31kpc (|Gillessen et al.ll2009D . determine 
the local circular speed v c (Rq) — RofiA* ~ Vq — (251 ±12 — 
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Figure 14. Full curves: the mean value of the thin-disc DF (upper 
panel) and the thick-disc DF (lower panel) over an ellipsoid in 
velocity space on which the model DF (13 1 I t is constant, divided 
by that constant value. Dotted curves: the rms variation over 
an ellipsoid in the model DF divided by the mean value on the 
ellipsoid. The x axis gives the semi-major axis of the ellipsoid in 
multiples of ay. 



V^jkms- MFlvnn et al.l (|2006T ) estimate the absolute I-band 
luminosity of the Galaxy to be Mi = —22.3. At this absolute 
magnitude the rid ge-line of the I-band Tully-Fisher relation 
i|Dale et al.lll999T) gives a circular speed of only 190 kms ; 
251 kms" 1 lies 2.6a from the ridge line. Thus the likelihood 
of the Galaxy in the context of the Tully-Fisher relation is 
small but increases rapidly with Vq, and this fact provides 
further support for an upward revision of Vq. 



3.4 Distinguishing the thin and thick discs 

Studies of the chemistry of the thick disc depend heavily on 
identifying nearby, bright stars that belong to th e thick disc 

llFuhrmann 



as ta r gets of medium-dis p ersion spec t rosco py 
19981; iBensbv et all 120031; IVenn et ail |2004 [I 



Bensbv et al 

20051 : iGilli et al.ll2006l ; iReddv et alj|2006t ). A popular strat- 
egy for identifying target stars is to assume ellipsoidal ve- 
locity distributions for each disc of the form l|Bensbv et al.l 
120031 ) 



(V - Va) 2 

2a 2 v 



W 2 



(31) 



where (U, V, W) are velocity components with respect to 
the LSR, and each component is assigned assumed values 
of au , av , aw and the rotational lag V a . A given star is as- 
signed to the population for which it gives the largest value 
of the DF that follows from equation (|31[) and assumed frac- 
tions of local stars that belong to each population. 
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Figure 15. Distributions in the (vr,v^) plane for v z = at three distances from the plane: 0.5 kpc (left), 1.5 kpc (centre) and 2.5 kpc. 
Exactly these distributions are predicted at z = and v z = 30, 67 and 97kms — 1 . 



The idea behind equation (|31[) is that, for appropriate 
parameters, / provides a useful approximation to the dfs 
of the two discs. Here we investigate the quality of this ap- 
proximation by comparing / with the thin- and thick-disc 
components of the sta ndard df when the parameters in / 
take the values used bv lBensbv et al.l(|2003D . which are given 
in Table [1] Fig. [14] quantifies the quality of the approxima- 
tion provided by / by reporting the fractional RMS variation 
in the appropriate component of the standard df over each 
ellipsoidal surface in velocity space on which / is constant 
(dotted curves) - ideally this would vanish. The full curves 
show the mean value of the df over an ellipsoid of constant 
/, divided by the value of / on that ellipsoid - since the nor- 
malisation of / is arbitrary, the full lines in Fig. [T3] may be 
shifted up or down at will, but must be horizontal if / it to 
provide a useful approximation to the relevant component 
of the standard df. 

The upper panel of Fig. [14] is for the thin disc and the 
lower panel for the thick disc. Since the full curve in the 
upper panel is approximately horizontal, we conclude that 
/ decreases from small to large ellipsoids in the same way 
that our model thin-disc df does. However, from the fact 
that the dotted curve lies above 0.7 for a > ou, we conclude 
that the model df varies by of order itself over the larger 
velocity ellipsoids of the thin disc. Thus / is a useful but not 
very accurate approximation to the model df for the thin 
disc. 

The lower panel of Fig. Q3] implies that / provides a 
very poor approximation to our model thick-disc DF: not 
only does the full curve rise by more than a factor of 4 as a 
increases to 2au, but the essentially linear rise of the dot- 
ted curve implies that the model df varies strongly over 
ellipsoids of constant /. What prevents / providing a good 
approximation to the df is the continuous increase in the 
asymmetric drift as the numerical value of the df decreases. 
This increase is apparent in the downward motion of isoden- 
sity contours in Fig. [TT]and drives the fall in (v<p) with z in 
the bottom right panel of Fig. 1101 In particular, at z — the 
density of stars in velocity space peaks close to the circular 
speed in the thick disc as in the thin. This fact conflicts with 
the structure of equation (pJTj) . 



3.5 Local and in-situ samples 

The standard df predicts that the distribution of stars in 
the (vrjV,/,) plane at height z and v z — is identical to the 



(vr,v^,) distribution of stars at any other height z < z and 
velocity v' z = {2[& z (z) — <£> z (z')]} 1/2 . Moreover, a sample 
of stars observed at some distance z from the plane will be 
heavily weighted towards stars whose vertical motions have 
turning points there. So there should be a close correspon- 
dence between local stars with a particular value of v z and 
samples of stars at a given value of z. Since the SDSS and 
its successors provide photometric distances and proper mo- 
tions for millions of stars that lie > 1 kpc from the plane, 
and the GCS catalogue provides space velocities for over ten 
thousand nearby stars, this correspondence can be tested in 
some detail. Fig. 1151 shows some sample (vn,v,p) distribu- 
tions. 



4 CONCLUSIONS 

We have explored the ability of distribution functions to pro- 
vide models of the thin and thick discs of the Milky Way. 
Our dfs are analytic functions of the actions of orbits, which 
ensures that there is an intuitive relation between the ob- 
servable properties of the population a df describes and the 
functional form of the df, and a meaningful way to com- 
pare models that use different gravitational potentials. In 
this paper we have used expressions for the actions that are 
only approximate, and imply that a star's vertical motion 
is adiabatically invariant during the star's motion parallel 
to the plane. In a forthcoming paper (McMillan et al., in 
preparation) orbital tori will be used to eliminate this ap- 
proximation, and thus quantify its validity. 

We have shown that the vertical density profile and 
kinematics of the disc are accurately modelled by the ex- 
tremely simple df ([9]). However, we rejected this df because 
it is essential to be able to break the df for the thin disc 
down at least into contributions from stars of various ages, 
and ideally into contributions from ranges in both age and 
metallicity. That is, we must recognise that the Galaxy is 
built up of innumerable stellar populations of various ages 
and metallicities, and each population has its own df. In this 
paper we have only begun to explore the resulting complex- 
ity by ascribing a single df to the thick disc and modelling 
the thin disc as a superposition of dfs for stars of different 
ages. In reality both discs are chemically inhomogeneous and 
we should assign a distinct df to the st ars born at each time 
with each chemical composition (e.g. Schonrich fc Binnevl 
2009). Hence the dfs presented in this paper should be con- 
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sidered building blocks from which more elaborate dfs may 
be in due course constructed. 

Our most basic building block is a "pseudo-isothermal" 
population of stars. Fig. Q] shows that the vertical distri- 
butions of young stellar populations is well modelled by 
a pseudo-isothermal population. The density of a pseudo- 
isothermal population does not decline exponentially with 
z, but Figs. l4l and 1 101 show that, remarkably, the composite 
population produced by stochastic acceleration of stars does 
have an exponentially decreasing density profile. An excel- 
lent fit to the observed density profile of the entire disc is 
obtained when a pseudo-isothermal thick disc is added to 
the composite thin disc. The dispersion in v z of thin-disc 
stars increases from 17.4kms _1 in the plane to 33kms _1 
at 2.5 kpc, while that of the thick-disc stars increases from 
~ 35kms _1 in the plane to ~ 48kms~ 1 at 2.5 kpc. The 
thick disc contributes to the solar cylinder 24 per cent of 
the luminosity contributed by the thin disc, or 19.4 per cent 
of the total luminosity of the disc. 

Even though we are assuming that the dynamical cou- 
pling between motions in and perpendicular to the plane is 
weak, two features of our dfs lead to strong correlations be- 
tween distributions in vr and v z . One feature is the fact that 
random velocities must increase as one moves inwards, and 
the other is the simultaneous increases in or and <r z that 
are driven by stochastic acceleration of a coeval population. 
Comparison of Figs. 171 and 1101 show that, on account of this 
correlation, the distribution of local stars in the (vr,v^) 
plane is atypical of the stellar population of the whole so- 
lar cylinder in just such a way that our composite disc df 
can simultaneously provide reasonable matches to the very 
different shapes of the distributions of GCS stars in vr and 
v$ . The widths of the model distributions in vr and v$ are 
controlled by a single parameter, a r o- The shape of the vr 
distribution is predetermined by our choice of the dfs func- 
tional form. The value of the parameter provides limited 
control of the shape of the v$ distribution and we obtain 
the best fit to the observed distribution when this param- 
eter is chosen such that the disc's surface density declines 
roughly exponentially with scale length 2.5 kpc, which hap- 
pens to agree with the sca l e leng th inferred from near-IR 
star counts bv lRobin et"afl i|2003h . 

In principle the df of the thick disc should be tightly 
constrained by the dependence on z of the velocity disper- 
sions or and a z . These dependencies have recently been 
determined for SDSS stars by two groups. Unfortunately, 
their results seem to be incompatible and the reasons for 
the conflict are unknown. 

The standard model provides an excellent fit to the sem- 
inal work of Kuijken & Gilmore, perhaps because the grav- 
itational potential in which the df is evaluated was partly 
fitted to that work. Some of the difficulties encountered here 
with fitting newer data may arise from inaccuracy of the po- 
tential used. A worthwhile exercise would be to fit data from 
the GCS, RAVE and SDSS surveys to models that combined 
dfs of the type presented here with and a multi-parameter 
gravitational potential: by simultaneously fitting the param- 
eters in both the df and the potential, one should be able to 
obtain reasonable fits to the data, providing the latter have 
been purged of such evident inconsistencies as those seen in 
Fig- El Data from more than one survey would probably have 
to be used since SDSS stars are too faint to constrain the 



thin disc tightly, although the RAVE survey, which certainly 
probes the thick disc effectively, may include enough nearby 
stars to make the Hipparcos-based GCS survey obsolete. 

The model fit to the distribution of GCS stars is far 
from perfect. Some of the disagreement arises because, as 
is well known, the Galactic bar and spiral arms give rise to 
features ("star streams") in the local velocity distribution 
that are inconsistent with the Galaxy being axisymmetric 
and in a steady state, as our models assume. Our favoured 
model v$ distribution would fit the data better if the con- 
ventional value of the solar motion Vq were ~ 6kms _1 too 
low. Tentative support for such an incre ase in the V"q is 
provided by astrometry o f stellar masers (jReid et al. I l2009l ; 
iMcMillan fc Binned 2009). and any increase would also tend 
to bring the Galaxy more into line with the Tully-Fisher re- 
lation between v c and Mj for external galaxies. By system- 
atically perturbing the velocities of all solar-neighbourhood 
stars, spiral structure might lead to the classical approach 
to the determination of Vq yielding an underestimate. Fur- 
ther work is required to explore this possibility, and at this 
stage we would merely stress that the systematic error in 
Vq is much larger than the formal errors given by DB98 
and lAumer fc Binnevl i|2009h . 

In the models, the asymmetric drifts of both the thin 
and thick discs increase with height. A disc's asymmetric 
drift is largely controlled by its parameter Rd and in the 
standard model the asymmetric drift of the thin disc exceeds 
that of the thick disc above 1 kpc because we have adopted 
a slightly larger value of Rd for the thick disc than for the 
thin disc. 

A popular strategy for assigning solar-neighbourhood 
stars to the thin or thick disc is to find the values taken 
by each disc's model df at the star's location. The model 
dfs used are perfectly ellipsoidal but we show that such dfs 
provide poor approximations to the thick-disc component of 
the standard df, so a markedly cleaner separation of the two 
discs could be obtained by replacing the ellipsoidal dfs by 
the thin- and thick-disc components of the standard df. 

Although the observational material relating to the 
Galaxy has increased enormously in recent years, we have 
shown that much of the available data can be successfully 
modelled with a simple analytical df. In a couple of aspects 
the data are in mild conflict with the df, but it is at least 
as likely that the fault lies with the data as the df. In the 
coming decade the volume and quality of the observational 
material available will increase dramatically. We anticipate 
that comparisons between each new data set and an evolv- 
ing standard df will reveal successes and failures similar to 
those encountered here. The successes will confirm the value 
of the df as a summary of a large and inhomogeneous body 
of data, and the failures will lead to critical re-examination 
of both data and df. Sometimes the failure will arise from 
a defective calibration of the data or incorrect assumptions 
used in its reduction, and other times it will indicate that 
the df is too simplistic. Either way we will learn something 
new and interesting. 

In this paper the df's parameters have been fitted to 
the data by eye and no attempt has been made to quantify 
uncertainties in parameter values. Clearly such uncertainties 
are important, and they could be most securely established 
by carrying the df's predictions closer to the raw observa- 
tions than we have done. In future work probability dis- 
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tributions in colour-magnitude-proper-motion space, etc., 
should be predicted that can be compared with the actual 
star counts. 

Upcoming infrared surveys, such as the VHS with Vista 
and APOGEE, will probe the disc at remote locations. The 
predictions of the standard df for those locations will be 
presented shortly, after orbital tori have been introduced as 
the means to convert between Cartesian and angle-action 
variables. This upgrade will make obsolete the approxima- 
tion of adiabatically invariant vertical motions used here. 



ACKNOWLEDGEMENTS 

I thank Michael Aumer for providing the data shown in 
Fig. [12] and Zeljko Ivezic for providing the data plotted in 
Fig. |9] The members of the Oxford dynamics group con- 
tributed valuable comments on drafts of this work. It is a 
pleasure to acknowledge valuable conversations with Ralph 
Schonrich. 



REFERENCES 

Abazajian K., et al., 2009, ApJS, 182, 543-558 
Ascasibar Y., Binney J., 2005, MNRAS, 356, 872 
Aumer M., Binney J., 2009, MNRAS, in press (arXiv 
0905.2512) 

Bensby T., Feltzing S., Lundstrom I., 2003, A&A, 410, 527 
Bensby T., Feltzing S., Lundstrom L, Ilyin I., 2005, A&A, 
433, 185 

Binney J., 1987, in "The Galaxy", eds. G. Gilmore & R. 
Carswell (Dordrecht: Reidel) p. 399 

Binney J., Merrifield M., 1998, "Galactic Astronomy", 
Princton University Press, Princeton 

Binney J., Tremaine S., 2008, "Galactic Dynamics", 
Princeton University Press, Princeton (BT08) 

Bond N.A., et al., 2009, arXiv 0909.0013 

Carlberg R.G., Sellwood J.A., 1985, ApJ, 292, 79 

Dale D.A., Giovanelli R., Haynes M.P., Campusano L.E., 
Hardy E., 1999, AJ, 118 1489 

Dehnen W., 1999, AJ, 118, 1201 

Dehnen W., Binney J., 1998, MNRAS, 298, 387 

Edvardsson B., Andersen J., Gustafsson B., Lambert DL., 
Nissen P.E., Tomkin J., 1993, A&A, 275 101 

Epchtein N., Simon G., Borsenberger J., de Batz B., Tan- 
guy F., Begon S., Texier P., Derriere S., and the DENIS 
Consortium, 2005, "DENIS Catalogue third data release" 
|http : / /cdsweb . u-st r asbg . fr /d enis . html | 

Flynn C, Holmber J., Portinari L., Fuchs B., Jahreiss H., 
2006, MNRAS, 372, 1149 

Fuchs B. et al., 2009, AJ, 137, 4149 (14 authors) 

Fuhrmann K., 1998, A&A, 338, 161 

Gillessen S., Eisenhauer F., Trippe S., Alexander T., Genzel 
R., Martins F., Ott T., 2009, ApJ, 692, 1075 

Gilli G., Israelian G., Ecuvillon A., Santos N.C., Mayor M., 
2006, A&A, 449, 723 

Gilmore G., Reid N., 1983, MNRAS, 202, 1025 

Gilmore G., Wyse R.F.G., Kuijken K., 1989, AnnRAA, 27, 
555 

Holmberg J., Flyn C, 2000, MNRAS, 313, 209 



Holmberg J., Nordstrom B., Andersen J., 2007, A&A 475, 
519 

Ivezic Z., et al., 2008, ApJ, 684, 287 

Juric M. et al., 2008, ApJ, 673, 864 

Kuijken K, Gilmore G., 1989, MNRAS, 239, 605 

Kuijken K, Gilmore G., 1991, ApJ, 367, L9 

Malhotra S., 1995, ApJ, 448, 138 

McMillan P.J., Binney J., 2009, MNRAS, submitted 
arXiv0907.4685 

Nordstrom B., Mayor M., Andersen J., Holmberg J., Pont 
F., J0rgensen B.R., Olsen E.H., Udry S., Mowlavi N., 
2004, A&A, 418, 989 

Oiling R.P., Dehnen W., 2003, ApJ, 599,275 

Perryman M.A.C., 1997, "The Hipparcos and Tycho Cat- 
alogues", (Noordwijk: ESA Publications) 

Reddy B.E., Lambert D.L., Allende Prieto C, 2006, MN- 
RAS, 367, 1329 

Reid M.J., Brunthaler A., 2004, ApJ, 616, 872 

Reid M.J., et al. (14 authors), 2009, ApJ, 700, 137 

Robin, A.C., Reyle, C, Derriere, S. & Picaud, S., 2003, 
A&A, 409, 523 

Schonrich R., Binney J., 2009, MNRAS, in press 

Shu F.H., 1969, ApJ, 158, 505 

Skrutskie M.F., et al., 2006, AJ, 131, 1163 

Spitzer L., 1942, ApJ, 95, 329 

Spitzer L., Schwarzschild M., 1953, ApJ, 118, 106 

Steinmetz et al., 2006, AJ, 132, 1645 

van der Kruit P.C., Searle L., 1982, A&A, 110, 61 

van Leeuwen F., 2007, Hipparcos, the New Reduction of the 
Raw Data, Springer Dordrecht 

Venn K.A., Irwin M., Shetrone M.D., Tout C.A., Hill V., 
Tolstoy E., 2004, AJ, 128, 1177 

Yanny B., et al., 2009a, AJ, 137, 4377 

Yanny B., et al., 2009b. larX"iv:0902. 17811 

Zacharias N., Urban S.E., Zacharias M.I., Wycoff G.L., Hall 
D.M., Monet D.G., Rafferty T.J., 2004, AJ, 127, 3043 



